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Using the liquid instanton model continued into Minkowski space and modified for finite 
temperature, the energy momentum tensor and the surface tension for buble walls during 
the QCD phase transition are derived. The resulting surface tension of bubble walls is in 
agreement with lattice calculations. Application to bubble collisions is discussed. 
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I. INTRODUCTION 



At the time 10 -5 s the temperature (T) of the universe was approximately 300 MeV. During the time 
interval between 10~ 5 -10 -4 s the universe passed through the temperature T c ~ 150 MeV, the critical tem- 
perature for the chiral phase transition from the quark-gluon plasma (QGP) to the hadronic phase (HP), the 
quark- hadron phase transition (QHT). There have been many numerical studies of the QHT. Some lattice 
gauge calculations indicate the the transition is weakly first order Q , which would imply bubble formation 
during the QHT. At the present time lattice calculations cannot determine the order of the transition 0]. 
Recent numerical calculations in the MIT Bag model || find a first order transition with hadronic bubbles 
stable for nucleation at a scale of 1 fm. 

Assuming that the QHT is first-order there have been a number of studies of the bubble nucleation Q 
and of possible observational effects of such a phase transition. The latter include the possibility of large 
scale density perturbations g and the seeding of primary magnetic fields || that could lead to observational 
effects in the Cosmic Microwave Background Radiation(CMBR), or large scale galactic or extra-galactic 
structure. 

The essential property of the bubbles needed to study the nucleation is the surface tension 0,0], o. The 
most recent lattice calculations give a value of cr ~ 0.015T C 3 g, while earlier studies |J find a larger value. 
In our present work we wish to estimate a from the QCD Lagrangian, with the result that we obtain the 
energy momentum tensor that can be used in future work on bubble collisions. For the electroweak phase 
transition bubble collisions have been studied and interesting estimates of magnetic fields formed during the 
collisions have been made [|lO|,^TJ using an abelian Higgs model plus the QED Lagrangian jl2|. However, 
because of the complexity of nonpcrturbative QCD, a reliable model of the bubble walls formed in the QHT 
has not been formulated. It is the goal of the present work to formulate such a model. 

Since the main structure of the bubble walls must be gluonic in nature we use pure gluodynamics. Noting 
that the walls separate the hadronic from the quark/gluon plasma regions of the universe, it is natural to 
consider an instanton picture of the walls. The QCD instantons are classical solutions for the color field, 
which were derived in Euclidean space ^3j,H. They connect points in two vacuua with different winding 
numbers. In analogy with Coleman's model jlq ] of regions of true and false vacuua, we consider the instantons 
connecting vacuua with different winding numbers in the quark-gluon region and the hadronic region on the 
opposite sides of the bubble. The Euclidean instanton model cannot be used for the treatment of bubble 
nucleation and collisions (in fact the energy density and surface tension vanish in this picture) , but the model 
must be continued to Minkowski space (JT5) . In brief, our picture is a Minkowski space analytic continuation 
of the instanton model for the bubble walls separating the quark/gluon from the hadronic regions during 
the QHT. 
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It is well known that the instantons of the instanton liquid model Jl6| can provide the essential nonpertur- 
bative gluonic effects of QCD at the medium length scale of 0.3-0.5 fm, although instantons do not provide 
confinement. As we show below, this seems to be the crucial length scale needed for QCD bubble walls. In 
recent work an effective Lagrangian was used to estimate domain wall formation possibly associated with the 
QCD phase transition |0. Although the resulting domain wall is within the hadronic bubble it resembles 
the bubble wall that we obtain with the instanton model in the present work. That model has also been used 
to examine the possibility of magnetic walls formed during the QHT seeding gallactic magnetic structure 
fist . Such a magnetic wall has also been considered for cosmic microwave background radiation polarization 
correlations fj"9[ |. Also, there are recent numerical investigations of the T-dependence of the properties of 
the QCD instantons which we need in the present investigation. 

In Sec. II we review the instanton liquid model at T=0 and discuss the energy density and equation of 
motion in Minkowski space. In Sec. Ill we use recent lattice results for instantons at finite temperature 
as well as our own recent work on glueballs |22j to estimate the surface tension. The surface tension 
is derived at time t=0, and serves as an initial condition for bubble nucleation and collisions in Minkowski 
space. In Sec. IV we discuss the energy momentum tensor for gluonic QCD in Minkowski space and possible 
applications based on the instanton model for bubble collisions and cosmological observations. 



II. INSTANTON MODEL AT T=0, CONTINUATION TO MINKOWSKI SPACE, AND QCD 

CHIRQAL PHASE TRANSITION BUBBLES 



In this section we review the instanton model of QCD at T=0, and the extension of the model to treat 
the bubble walls in the QCD chiral phase transition. Instantons are obtained and studied in Euclidean 
space. As described in early work on the application of instantons to bubble nucleation and collisions |l5| ] 
in formulating the instanton bubble itself one proceeds in Eucledean space, while for the study of collisions 
of nucleating bubbles one must work in Minkowski space. This will be discussed further in Sec. IV. 



A. T=0 Instanton Model in Minkowski Space 



The Lagrangian density for pure glue is 



where 



£glue _ —-G ■ G, (1) 



G" u = d„A v - dyAp - ig[A^ A,} (2) 
A, = A$\»/2, 



with A™ the eight SU(3) Cell-Mann matrices. The instantons are the classical solutions of the gauge fields 
for a pure SU(2) gauge theory. I.e., w 
field, if one keeps only the instanton glu 
the classical action in Euclidean space 



for a pure SU(2) gauge theory. I.e., writing ^4™ = A™' lnst + A™' qu , where A™' mst is a pure gauge classical 
field, if one keeps only the instanton gluonic field , and solves the equation of motion obtained by minimizing 



6[ / G mst ■ G mst ] = 0, (3) 



one obtains the solution 11131 



(.r-2 _|_ p2\ 



Al' mst {x) = 



n "4o 2 
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for the instanton and a similar expression with -n for the anti-instanton, where p is the instanton size and 
the 77™,, are @: r\l c = e abc ;ri^ 4 = 5 a ^?]%„ = -5 au , with (a,b,c) = (1,2,3) and (n,v)= (1...4). The instanton 
connects points in two QCD vacua which differ by one unit of winding number. The instanton contribution 
to the Lagrangian density is 
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-G mst ■ G mst = 48 , , P OXA . (5) 

The energy momentum tensor, T^, can be obtained from the Lagrangian, and in Minkowski space is 
given by (see, e.g., Ref |§) 

= ^(GT^L - \g^GfG a0a ). (6) 

a 

To get the Euclidean space expression for the energy density, T 44 , one carries out the analytic continuation 
giving x 2 = x^x^ — x 2 + x\ in Euclidean space, while it is x 2 — Xq in Minkowski space. From this one see 
that g 44 = — g 00 . Therefore in Euclidean space the energy density in the instanton model is 

Tiirtst f 1 s~i 1 ~^ /~i r ' 

44 — " (1 4q(j (1q 4 -f — <-T act/ 3Lz aQ ^ 

= 0. (7) 

The result shown in Eq. (^) follows from Eq. (||) and the properties of the 77™^ given in Ref [jl4j . 

For the nucleation of QCD bubles and for bubble collisions one must work in Minkowski space. From 
Eqs.(||4|) one finds for the energy density in Minkowski space 

rpQ0,i7ist ^_Qinst Qinst 

~ 2 

= 96 7 2u ( g ) 
(x z + p z p 

Eq.(||) gives the energy density for an instanton. To complete the instanton model for T 00 one must 
find the instanton density, n, for the particular system in question. Evaluation of the instanton density is 
reviewed for both T=0 and finite T in Ref ]l6[| . In the instanton model the density is easily evaluated in 
terms of the gluon condensate, which is proportional to the vacuum matrix element < G ■ G >. Thus the 
instanton vlaue for n, n(inst), is given by 

/ - \ <i G • G ^> , . 

n{inst) = — - — — — -. (9) 

From Eq.(|j) and the phenonomenological value for the quark constant < G ■ G >~ l.OGeV 4 , one finds that 

n(inst) ~ l.O/VnT 4 . (10) 

This value is consistent with that determined by the tunneling amplitude ||l4[ , |l6|| . Another evaluation of 
n can be obtained from the QCD sum rule evaluation of scalar glueballs. Since the correlator for scalar 
glueballs is given by the Fourier transform of < G(x) ■ G(x)G(0) ■ G(0) >, the pole term in the dispersion 
relation for the correlator has a factor of n. By comparing to the 1.5 GeV candidate for a scalar glueball, 
the best fits were found 0] for values up to about twice the instanton value. Also, some lattice calculations 
|20f have found larger values. 

To evaluate the surface tension of the bubble walls we must also consider instantons at finite temperature 
during bubble nucleation, which we do in the next section. 
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B. QCD Chiral Phase Transition Bubble Walls 



The equations of motion for pure glue QCD are obtained from the Lagrangian given by Eq. (jl|)- Using 

5[[g-G] = 0, (11) 



and working in a SU(2) version of QCD, with — A^a a /2, er° being the Paili spin matrices, and [A^, A u ] = 
iCabcA^A^a 11 /2, the equations of motion are 

d^Al + ge ahc (2A b fl d fl A c v - A\d v A^ c ) + g 2 e abc e cei ' A\A» e A[ = 0. (12) 

Recognizing that the instanton form of Eq. (^) is a classical solution to these equations of motion, we try the 
form 

A n / nst (x) = 2 V -^F(x 2 ), (13) 
and sustituting in Eq.(H2h find that F(x 2 ) satisfies 



dF 

d 2 F = -4— - 12F 2 + 8x 2 F 3 (14) 
ox z 

in the Lorentz gauge, with the gauge condition d^A^ — 0. By using the initial condition 

F(x,t = 0)= J- 2 , (15) 
[x - Xo) + P 

one finds that the instanton form is indeed an approximate solution. This is discussed further in Section IV, 
where collisions of QCD bubbles are discussed. 

From the phenomenological applications [ fl6| the instanton liquid model was developed, with p = 0.33 fm. 
This gives the thickness of the instanton wall at temperature T=0. To complete our model we need the 
instanton density at T=T C , as well s the value of the instanton size at that temperature. This is discussed 
in the next section. 



III. BUBBLE WALL SURFACE TENSION 



In this section we evaluate the surface tension in the instanton model. In order to do this we must consider 
the dependence of the instanton size and density on the temperature. We start with calculating the surface 
tension for an instanton wall. 

The surface tension is the essential parameter for classical bubble nucleation. For the treatment of nucle- 
ation one must work in Minkowski space. As discussed in detail in Ref. |l5[ , to treat nucleation and collisions 
with instanton input one makes an analytic continuation from Euclidean to Minkowski space at the initial 
time, t=0. For our problem this is achieved by replacing t by it, or in Minkowski space x^x^ — x 2 — t 2 . In 
the present section we study the surface tension at t=0; and the energy density in Minikowski space derived 
in the previous section from the instanton solutions can be used for the derivation of the surface tension. 
The results at t=0 which we find will serve as part of the initial conditions for the study of nucleation and 
collisions, as discussed in Sec. IV. 

The surface energy is the surface tension x the area of the wall for a thin wall. Note that the wall 
thickness is of the order of p, the instanton size, which is a fraction of a fermi, so that a thin-wall formalism 
is appropriate. Therefore the surface tension with one instanton at t=0 is 
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= /^°°*"M,o) 



From Eqs.dQ 



inst = 3 . 2 o I dx 



dx-G-G. (16) 



5 / J_ 



(X 2 +p 2 ) 4 



!?• < 17 > 

In order to complete our model, we must work with instantons at the temperature T = T c , which means 
using a modified instanton size and density. Moreover, since we are dealing with a gluonic wall separating the 
two phases, the QGP and the HP, the number density is quite different from the result given in Eqs.(p[lC|) 
for the T=0 system. We use the model applied for the lattice calculations ||, and described in detail in 
Rcf. ||. In this picture the free energy is calculated using two peaks for the two phases on the two sides 
of the bubble, separated by a mixed quark/gluon-hadronic phase. There will be a factor of the instanton 
number density, i.e., the instanton density times the effective four- volume, for each peak. Therefore, from 
from Eq.()l7j) the resulting surface tension in our instanton model is 

a = 6-5 (18) 
V 

where N is the instanton number density at the surface and p is the instanton size at T = T c . We obtain N 
from the T = instanton number per four-volume, n = (N/V) of Eq.(]9j). Defining n = n GeV 4 at finite T 
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and using V= 7r\ N = 5.96 n 2 , where we use p = 0.25 fm, as discussed next. 

At T=0 p = 0.33 fm in the liquid instanton model j|6). This has also been found in lattice calculations 
p4| . There have been a number of numerical investigations on the modification of the instanton system at 
finite T. Based on these calculation [^0|,^l) we take p = 0.25 fm. There is uncertainty in the value of the 
instanton density, n, at T=0. In our work on scalar mesons/glueballs |^2| we found values of n = 0.0008- 
0.0015, with 0.0008 used in the instanton liquid model. These values also cover the range of values of lattice 
gauge calculaitons over the past decade. From the results of the finite T calculations the n seems to decrease 
by about a factor of 2 at T c , so we use the range n — 0.0004-0.00075. This gives us the result 

a = (0.013 -> 0.046)T C 3 , (19) 

compared to the value a = 0.015 T 3 in the most recent lattice calculation ||. We note that the instanton 
liquid model of the QCD bubble wall gives a value a = 0.013 Tj? in very good agreement with the lattice 
calculation, and that our model for the bubble wall is quite promising. 



IV. DISCUSSION OF APPLICATIONS TO COLLISIONS 



Although the knowledge of the surface tension is essential for studing the nucleation and collisions of 
bubbles during the QCD chiral phase transition, there is no observable directly related to these bubbles, so 
the challenge is to find observable effects related to these bubbles. In our recent work we have investigated 
the possibility of a magnetic walls being formed by interactions of hadrons with an internal gluon instanton 
wall which might be formed within the hadronic phase during the QCD chiral phase transition. In this 
section we discuss this possibility and the relationship to the wall surface tension that has been derived 
above. 
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As was discussed in Sec. II the instanton model itself describes a static bubble. For the study of bubble 
nucleation and collisions one works in Minkowski space. Starting from the instanton Lagrangian and using 
standard field theory methods one obtains the color E and B fields, 



77m /~in 

&i — u i0 



B? = 



\e ijk G% (20) 



with the i,j,k indices running 1,2,3. From this the energy-momentum tensor, T for an instanton in Minkowski 
space is given as [E3| 



nOQ,inst 



1 (E™ ■ E™ + B" ■ B") (21) 
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rpOi,inst c^^E n B n (22) 

rpij.inst (E nl E n ^ -(- B nl B n ^ rp00,inst g (23) 

From the results of Sec. II and with a straightforward calculation we obtain for the stress energy tensor in 
the instanton model from Eq.(|23|) 



2 

(x 2 + p 2 ) 2 , 



T i jM st = _ ( . ^P 2 _ 2 ^ ) s ti . (24) 



To obtain the color E abd B fields one must solve the QCd equations of motion. Solutions to the equations 
of motion with an effective Lagrangian have been found for a pure SU(2) gauge theory in Minkowski space 
pfjj. Using the general picture of Coleman Jl5| in which one sets the instanton initial conditions at time 
t=0 and then uses the equations of motion in Minkowski space for the evolution, a SU(2) model of of QCD 
instantons has been developed and applied to heavy relativistic heavy ion collisions I27J . 

We are developing a similar program for our cosmological studies of the QCD chiral phase transition in 
the early universe based on QCD field theory. During the collision of two bubbles the colliding walls might 
form an internal wall within the merged bubbles. In classical nucleation theory when two identical bubbles 
collide they merge to form a bubble with an interior wall having the same surface tension as the original 
bubble walls. In our case the result would be instanton walls described in the previous sections. Note that 
in Ref. |l7j an estimate of the lifetime of the gluonic domain walls in their model was 10~ 5 s, which allows 
magnetic walls to be formed Although the present theory is quite different, a two-dimensional model 
of bubble collisions starting from the gluonic QCD Lagrangian finds that an internal gluonic wall is formed 
pifl , and this suggests that magnetic structures could form. Using the standard evolution methods as used 
for tangled magnetic fields for the scattering (25) or for metric perturbations [^9| one can investigate possible 
observational effects in the CMBR. A preliminary investigation of this has been carried out assuming the 
interior instanton wall is formed, which is shown to lead to an interior magnetic wall, and CMBR polarization 
correlations are found that could be observed in the next few years. 

The crucial question is whether an internal wall is formed with the properties of the nucleating QCD 
bubbles. In a 1+1 dimensional study in pure QCD, the equations of motion of Eq.(|l4|) were solved. This 
corresponds physically to very large bubbles colliding, as illustrated in Fig.l. 
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Fig. 1 <^>C M > bubble walls colliding 

To model collisions one takes an initial condition with two walls, such as 

F(x, t = 0) = -j + 1 2 , (25) 

with 2xo the separation between the two bubble walls at time t=0, and follows the evolution. In Ref ps| ] 
it was found that an interior wall is indeed formed, which seems to be similar to the initial walls. Fig.l 
illustrates the process at the instant when the walls are overlapping, and the calculation shows the two 
bubbles separating with an internal instanton-likc structure remaining at the collision region. From this one 
would conclude that the interior wall would have approximately the same surface tension as discussed in the 
previous sections, and could result in the magnetic wall of Ref Jl9| . Moreover, the studies of the QCD chiral 
phase transition find that in contrast to the electroweak phase transition where many bubles nucleate, 
the QCD chiral phase transition seems to proceed via inhomogeneous nucleation, with larger distance between 
bubbles and few nucleatingbubbles involved in the transition to the hadronic phase. Thus the large bubble 
collision of Fig.l and Ref J2S| is well-founded. 

To investigate the nucleation and collision of bubbles during the QCD chiral phase transition one must 
include quark and hadronic degrees of freedom, since they play a major role in the free energy difference 
between the hadronic and quark/gluon phases. With Eq. (|i~2]) , however, one can learn a great deal about the 
bubble walls, and our results for the surface tension in the present paper are quite promising. The full 3+1 
dimensional treatment of nucleating QCD bubbles is a subject of future research. 
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